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SECTION  I 


INTRODUCTION 


Several  questions  have  been  raised  regarding  the  propagation  of  air 
blast  through  tunnels  with  various  input  and  surface  conditions.  A 
series  of  experiments  was  conducted  to  investigate  this  problem  during 
the  Pre  Dice  Throw  100  ton  TOT  event  by  TRW  Systems.  In  addition  TRW 
and  the  Air  Force  Weapons  Laboratory  (AFWL)  have  made  several  hydro- 
dynamic  calculations  in  support  of  the  experiments  and  the  underlying 
question  of  blast  propagation  in  general. 

The  subsequent  sections  of  this  report  include  a summary  of  the 
AFWL  calculations. 


SECTION  II 


SMOOTH  TUNNEL  PROPAGATION 

Four  calculations  have  been  completed  for  a smooth  walled  tunnel.  In 
each  case  the  input  boundary  conditions  at  the  tunnel  entrance  were  taken 
from  the  AFWL  1 KT  Standard  (ref.  1)  scaled  to  1 MT.  The  calculations 
were  carried  to  a distance  of  3.2  km,  which  corresponds  to  a free  air 
overpressure  of  less  than  10  psi.  The  inlet  peak  overpressure  levels 
were  10  , 10  , 600,  and  100  psi.  Table  1 contains  the  initial  conditions 
for  each  of  the  four  calculations.  All  calculations  were  made  with  a 
one-dimensional,  Eulerian,  second  order  accurate,  hydrodynamics  cO'de 
(DYTSAP) . The.  zone  size  in  the  vicinity  of  the  shock  front  was  held 
constant  at  2 meters.  Figures  1 through  4 are  plots  of  the  peak  over- 
pressure versus  distance  from  burst  point  for  the  four  inlet  pressure 
levels  above. 

The  1 KT  standard  used  as  boundary  conditions  is  an  analytic  fit 
with  discontinuous  shocks.  Hydrodynamic  codes  define  the  shock  over  a 
minimum  of  three  zones,  thus  the  shock  front  for  these  calculations  is 
spread  over  approximately  6 meters  (m) . The  effect  of  this  interface 
between  analytic  and  zone  definition  is  an  artificial  drop  in  the  peak 
overpressure  immediately  upon  entering  the  hydro  grid.  As  a check  on 
this  phenomenon,  an  additional  calculation  was  made  with  0.06  m zones 
(figure  5) . The  fine  zone  calculation  exhibited  the  same  drop  in 

1.  Needham,  C.  E. , Air  Force  Weapons  Laboratory  (1  KT)  Nuclear  Blast 
Sjtarrdard^,  AFWL-TR- 73-55,  Air  Force  Weapons  Laboratory,  Kirtland  AFB 
NM,  April  1975. 
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Table  1 


Problem 

Yield 

Starting  time  (sec) 

Distance  of  left 
mesh  edge  from 
burst  center  (H) 

Initial  peak 
pressure  (Pascals) 


Number  of  cells 
Width  of  cells  (M) 


1-D  INPUT  PARAMETERS 
13.0011  13.0012 


5.1  x 10' 


- 1 


9.3  x 10 


-2 


9.39  x 102 

6.8  x 105 
(Approx. 
100  psi) 

900 

2. 


4.54  x 102 

4.129  x 106 
(Approx. 

600  psi) 

900 

2. 


13.0013 

1 

7.3  x 10"3 

1.694  x 102 

6.8  x 10s 
(Approx. 

10*  psi) 

900 

2. 


13.0014 

1 

1.0  x 10“3 

7.96  x 101 

6.5  x 10® 
(105  psi) 

900 

2. 
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OVERPRESSURE  (PASCALS) 


11 


E (PASCALS) 


DISTANCE  ( METERS) 


DISTANCE  (METERS) 


peak  overpressure  but  only  a small  percentage  of  that  of  the  2 m zone 
calculation.  As  the  two  calculations  progress,  the  overpressure 
distance  curves  merge,  and  by  a distance  of  500  m no  significant 
differences  can  be  found.  The  conclusion  is  that  the  correct  amount 
of  energy  is  entering  the  calculational  mesh  but  takes  a small  amount 
of  time  (about  20  zones)  to  adjust  to  the  new  conditions. 

Table  2 compares  the  various  results  with  one  another  and  with  the 
free  air  case  at  1 and  3 kilometers  from  the  burst  point. 
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Table  2 


OVERPRESSURE 

FOR  VARIOUS 

INLET  CONDITIONS 

P 

(Pascals) 

P (Pascals) 
at  R - 1 km 

P (Pascals) 
at  R = 3 km 

6.8  x 108 
(10s  psi) 

1.2  x 107 

3.5  x 106 

6.8  x 10 7 
(104  psi) 

4.3  x 106 

1.2  x 106 

4.13  x 10s 
(600  psi) 

1.0  x 106 

3.2  x 10s 

6.5  x 10s 
(100  psi) 

4.9  x 10s 

1.5  x 10s 

Free  air 

4.8  x 10s 

6.0  x 104 

16 


The  table  shows  that  a factor  of  three  increase  in  pressure  results  at  a 
distance  of  3 kilometers  when  the  inlet  pressure  is  increased  from  104 
to  10s  psi.  The  pressure  at  3 km  for  104  inlet  pressure  is  20  times  as 
great  as  that  of  free  air  at  the  same  point.  The  overpressure  impulse 
is  more  than  200  times  that  of  free  air.  Similar  comparisons  can  be 
made  using  table  1 for  other  inlet  pressures.  Plots  of  the  pressure 
waveforms  for  these  four  problems  are  contained  in  Appendix  A.  A tabula- 
tion of  inlet  parameters  for  the  104  and  10s  psi  waveforms  is  given  in 
Appendix  B. 
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SECTION  III 


THE  TWO-DIMENSIONAL  CALCULATION 

One  calculation  was  made  in  two-dimensional  cylindrical  geometry  to 
study  the  effects  on  air  blast  of  corrugations  on  the  inner  surface  of 
the  tunnel  wall.  The  calculational  mesh  was  28  radial  by  800  axial 
zones  for  a total  of  22,400.  The  radial  zone  dimensions  were  6 cm  near 
the  outer  edge  of  the  tunnel  with  gradually  increasing  size  to  15  cm  at 
the  center  of  the  tunnel.  The  800  axial  zones  were  initially  12  cm,  and 
the  problem  wras  rezoned  such  that  the  shock  front  always  traveled  through 
12  cm  zones.  The  tunnel  had  an  outside  diameter  of  4.68  m.  Attached  to 
the  wall  were  perfectly  reflecting  rings  extending  0.24  m into  the  tunnel 
and  0.6  m of  outer  wall  between  rings  or  2.04  m from  center  to  center. 

At  the  tunnel  entrance  lO4  psi  input  from  a 2-MT  free  air  burst  was  used 
as  a boundary  condition.  The  calculation  has  been  run  a distance  of  180+ 
meters  down  the  tunnel. 

Waveform  plots  indicate  multiple  shock  enhancement  on  the  axis  of 
symmetry.  Overall  decay  of  the  shock  front  pressure  appears  to  be 
matching  that  of  the  free  air  curve  and  is  significantly  more  rapid  than 
the  1-D  (smooth  walled)  calculation.  The  enhancement  of  reflected  shocks 
is  a significant  factor.  A train  of  shock  waves  of  gradually  decreasing 
peak  pressure  follow  the  main  shock  at  a nearly  constant  time  interval 
between  1 and  1.5  ms.  The  impulse  is  at  least  as  great  as  that  of  the 
1-D  calculations. 
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The  inlet  temperature,  after  shock  passage,  is  in  excess  of  25,000°K. 

This  high  temperature  gas  is  pulled  into  the  tunnel  behind  the  shock  and 
is  further  shock  heated  to  temperatures  over  30,000°K. 

Table  3 gives  the  maximum  pressure  reached  at  several  points  in  the 
tunnel.  In  almost  all  cases  the  peak  pressure  reached  is  not  the  first 
shock  but  the  second  or  third.  Figure  6 is  the  pressure-distance  curve 
inside  the  corrugated  tunnel.  Also  included  are  curves  from  the  free 
air  case  and  the  smooth  tunnel  calculation.  Appendix  C contains  a 
number  of  plots  of  a small  portion  of  the  length  of  the  tunnel  during 
shock  passage.  The  multiple  shock  enhancement  is  clearly  shown. 

1 

8 

I 

v 
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Table  3 


PEAK  PRESSURES  IN  CORRUGATED  TUBE 
^angG  Peak  pressure  (dynes /cm2 I 


burst 

(meter) 

Tunnel 

center 

Between  center 
and  outer  edpe 

Near  out- 
side edge 

214.3 

6.77E8 

6.80E8 

1.05E9 

224 

7.81E8 

3.85E8 

4.43E8 

234 

5.55E8 

3.80E8 

3.97E8 

244 

4.40E8 

3.85E8 

3.61E8 

253 

3.76E8 

3.32E8 

4.74E8 

263 

3.15E8 

2.96E8 

3.33E8 

273 

3.21E8 

2.66E8 

2.41E8 

284 

3.17E8 

2.60E8 

6.04E8 

294 

2.89E8 

2.32E8 

3.71E8 

304 

2.43E8 

2.00E8 

2.43E8 

314 

2.84E8 

2.53E8 

2.64E8 

324 

2.54E8 

2.14E8 

1.78E8 

334 

2.13E8 

1.76E8 

1.22E8 

344 

1.83E8 

1.47E8 

1.17E8 

353 

1.40E8 

1.36E8 

1.76E8 

363 

1.31E8 

1.25E8 

5.87E7 

RANGE  (METERS) 


APPENDIX  A 


I-D  WAVEFORM  COMPARISONS 


This  appendix  contains  comparisons  of 
the  tunnel  with  free  air  waveforms  at  the 


overpressure  waveforms  inside 
corresponding  times. 
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APPENDIX  B 


TABULATED  INLET  CONDITIONS 


This  appendix  contains  tabulated  values  of  the  inlet  conditions  for 

4 5 

the  10  and  10  psi  conditions.  All  units  are  CGS . 
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APPENDIX  C 

PLOTS  OF  THE  2-D  CALCULATION 

This  appendix  contains  histograms,  vectors,  and  contour  plots  of 
pertinent  variables  foi  the  corrugated  tunnel  calculation.  The  shocks 
reflected  from  the  corrugations,  one  shock  for  each  corrugation,  are 
clearly  shown.  The  vortex  formed  between  the  corrugations  can  also  be 
seen  in  the  velocity  vector  plots  and  on  the  vorticity  plots.  The 
multiple  shock  fronts,  so  clearly  seen  on  the  pressure  histograms,  are 
not  numerical  instabilities.  This  is  demonstrated  by  the  fact  that  in 
the  first  eight  shocks  nearly  100  zones  are  used  to  resolve  them, 
instabilities  would  be  seen  in  essentially  every  other  zone. 

I 

* 

The  tunnel  is  oriented  along  the  axis  of  symmetry  of  the  code. 

Thus,  on  the  plots,  altitude  refers  to  distance  along  the  length  of 
the  tunnel  from  the  opening,  not  from  the  burst  center.  The  numbers 
on  the  right  and  the  top  of  the  graphs  are  the  zone  numbers  for  the 
calculational  grid. 
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221  N.  Grand,  St  Louis,  MO 
63100 
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